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ABSTRACT 

We measure apparent velocities {v app ) of absorption lines for 36 white dwarfs (WDs) with helium- 
dominated atmospheres - 16 DBAs and 20 DBs - using optical spectra taken for the European 
Southern Observatory SN la progenitor survey (SPY). We find a difference of 6.9 ± 6.9km s _1 in the 
average apparent velocity of the Ho lines versus that of the He I 5876 A for our DBAs. This is a measure 
of the blueshift of this He line due to pre ssure effects. By us ing this as a correction, we extend the 
gravitational redshift method employed bv lFalcon et al.l (|2010f ) to use the apparent velocity of the He I 
5876 A line and conduct the first gravitational redshift investigation of a group of WDs without visible 
hydrogen lines. We use biweight estimators to find an average apparent velocity, (u a pp)Bi, (and hence 
average gravitational redshift, (v s )bi) for our WDs; from that we derive an average mass, (M )bi- For 
the DBAs, we find (i> a pp) B i = 40.8 ± 4.7km s" 1 and derive (M) B i = O.7lt o ;^M . Though different 
from (f a pp) of DAs (32.57 km s _1 ) at the 91% confidence level and suggestive of a larger DBA mean 
mass than that for normal DAs derived using the same method (0. 647+^014 M^ lFalcon et al.ll201Q[ ). we 
do not claim this as a stringent detection. Rather, we emphasize that the difference between (i> app )Bi 
of the DBAs and {v app ) of normal DAs is no larger than 9.2 km s _1 , at the 95% confidence level; this 
corresponds to roughly 0.10 Mq. For the DBs, we find (i>^ p )bi = 42.9 ± 8.49 km s" 1 after applying 
the blueshift correction and determine (M}bi = 0.74toog^o- The difference between (w a i p c p )Bi of the 
DBs and {v app ) of DAs is < 11.5km s _1 (~ 0.12 M@), at the 95% confidence level. The gravitational 
redshift m ethod indicates m uch larger mean masses than the spectroscopic determinations of the same 
sample bv lVoss et al.1 (|2007t ). Given the small sample sizes, it is possible that systematic uncertainties 
are skewing our results due to the potential of kinematic substructures that may not average out. We 
estimate this to be unlikely, but a larger sample size is necessary to rule out these systematics. 
Subject headings: stars: kinematics and dynamics - techniques: radial velocities - techniques: spec- 
troscopic - white dwarfs 



1. INTRODUCTION 



In lFalcon"etlal([2lHll . we show that the gravitational 
redshift method is an effective tool for measuring mean 
masses of groups of white dwarfs (WDs) and has the ad- 
vantage of being m ostly independent from the spectro- 
scopic method (e.g.. iBergeron et aLlll992f ). iFalcon et al.l 
(|2010f) investigate normal DAs, the largest class of WDs. 
The next logical step is to ask whether the method can 
be applied to the second largest class, DBs, which con- 
stitute 20% of all WDs below T n ff ~ 17,000K and ~ 9% 
of WDs at higher tem peratures (jBeauchamp et al1ll996t 
IBergeron et~aljl2011| ). 

However, the gravitational redshift method histori- 
cally uses the apparent velocities of hydro g en Ba lmer 
lin e cores. The wor k by | Shipman k, Mehanl (|1976l ) and 
by iGrabowski et al.l (|1987f ) show Ha to be suitable for 
this purpose since it are not significantly affected by 
pressure shif ts. Pure DB spectra e xhibi t only helium 
lines, and as iGreenstein &; Trimble! ()1967l ) first pointed 



out, using WD photospheric helium lines for gravita- 
tional redshift measurements can be difficult due to 
the likelihood of systematics introduced by pressure ef- 
fects. These systematics include tha t, in theory, and 
with some experi mental support (e.g.. iBerg et al.lll962t 
I Perez et al.l I2003D , different helium lines can be pres- 
sure shifted by different amounts, in different (blue or 
red) direct ions, and w i th a dependency on tempera- 
ture (e.g.. iGriem et al.1 11962) : iBassalo fc Cattani 119761 : 
iDimitrijevic fc Sahal-Brechoti 119901 : lOmar et al.l l2006[ ). 
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For this reason, attempts at gravitational redshift mea- 
surements for he lium-dominated WDs have been sparse. 

iKoesterl (|1987t) measures line shifts of He I 4026, 4471, 
4713, and 4922 A in the spectrum of the common proper 
motion star WD 0615-591 and of He I 4471 and 4922 A 
in the wide binary WD 2129+000. These line shifts 
are negative (blue), meaning that the magnitude of the 
pressure effects are larger than the magnitude of the 
gravitational redshift, which, in this case, are opposing 
each other. The fact that these WDs are relatively cool 
(|Bergeron et al.|[201 If ) is consistent with the expectation 
that pressure effects should be significant; we will elab- 
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orate on this point in Section 13.61 iKoesterl (|1987[ ) con- 
cludes that due to the state of the theory at the time - 
laboratory measurements and theoretical predictions of- 
ten disagreed on the magnitude and sometimes even sign 
of the shift - he cannot deduce meaningful gravitational 
r edshifts for these two DBs. 

IWegner et all £[989) measure the gravitational redshift 
for the Hyades DBA WD 0437+138 using Ha and men- 
tion that the velocity "...is unaffected by the pressure 
shift problems of helium" while providing no further 
detail. The importance of van der W aals broadening, 
in co ol helium-dominated atmospheres (jBergeron et al.l 
119911) . however, was perhaps not yet well-established. In 
hindsight, it is likely that Ha in this WD is significantly 
affected. 

The sample of common pr oper motion binary and 
cluster WDs in Rcidj (]1996f ) contains three targets 
with helium-dominated atmospheres. For both DBAs, 
WD 0437+138 and WD 1425+540, Reid determines dif- 
ferent gravitational redshifts from using Ha than from 
H/3. He mentions that this discrepancy could be because 
of pressure shifts due to the high atmospheric helium 
abundance and deems the Ha result as the better red- 
shift estimate since this line sh ould be l ess aff ected than 
H/3. Reid also measures, like IKoesterl l(l98l . negative 
shifts of helium lines for the DB WD 2129+000. One of 
his measured lines ( He I 4921, 5015 , and 6678 A) is in 
common with that of lKoesterl (|1987D . 

With more recent high-resolution sp ectroscopic obser- 
vations of helium- atmosp here WDs (Voss et al.l 120071) 
and with the method of iFalcon et al. ( 2010D . we now 
have the tools to revisit the gravitational redshift of 
DBs. Such an investigat ion is a valuable check to the lat- 
est spectroscopic work (jBergeron et al.l f201 lh . the anal- 
ysis of which is nontrivial due in part to the chal- 
lenge of interpreting p r essure-broadened helium lines 
(IBeauchamp et alj|1997t iBeauchamp fc Wesemaell 119981: 
iBeauchamp et al.lll999[ ). 

In Section 13.61 we discuss using the apparent veloc- 
ity of t he He I 5876 A l ine in the context of our work 
and, as IWegner fc Reidl (|1987l ) suggest, check for consis- 
tency within the DBA sample of the apparent velocities 
of both the hydrogen and helium line species. In Section 
I4.1l we perform the original gravitational redshift method 
that uses hydrogen Balmer lines on the DBAs, which re- 
sults in the most d irect comparison with the DAs from 
IFalcon et all (|2010D . Then in Section [O] we extend the 
method to WDs with only helium lines. 

2. OBSERVATIONS 

We use spectroscopic data from the European South- 
ern Observatory ( ESQ) SN la progenitor survey (SPY; 
iNapiwotzki et all [2001). These observations, taken 
using the UV- Vi sual Echelle Spectrograph (UVES; 
iDekker et al.l 1200 (1) at Kueyen, Unit Telescope 2 of the 
ESO VLT array, constitute the largest, homogeneous, 
high-resolution (0.36 A or ~ 16 km s _1 at Ha) spectro- 
scopic dataset for WDs. We obtain the pipeline-reduced 
data online through the publicly available ESO Science 
Archive Facility. 

2.1. Samples 

The helium-atmosphere WDs in our sample are derive d 
from the SPY objects analyzed by IVoss et all (|2007l ). 



Out of their 38 DBAs with atmospheric parameter deter- 
minations, we exclude one magnetic WD, one DO, and 
twelve coo l WDs with T e s < 16,500 K. For these twelve, 
IVoss et al~l (j2007| ) fix the surface gravity at log g=8 in 
order to obtain spectral fits; this forfeits our ability to 
properly compare these stars with our gravitational red- 
shift results. Due to data quality, we are unable to satis- 
factorily measure apparent velocities (w a pp) of Ha for an 
additional eight, which leaves us with 16 DBAs. Keep in 
mind that for these objects, the atmospheric abundance 
of hydrogen is small and Ha is often barely visib le. 

Out of the 31 measured DBs from IVoss et al.1 (|2007| ). 
we exclude ni ne cool WDs and two known to have cool 
companions (iZuckerman fc Becklinl 119921: iFarihi et all 
120051: iHoard et al.l 120071 ). We end up with 20 normal 
DBs and 36 helium-atmosphere WDs overall. 

Four of our DBAs show Ca II lines in their spec- 
tra, giving them the additional classifi cation as DBAZ 
([Koester et al.l l2005bl: IVoss et all 120071) and two DBs 
are classified as DBZs ( Voss et al.l 120071 ). Since the 
SPY Collaboration has not reported any of our WDs 
to have stellar companions, we presume them to not 
be in close binary systems, though two of our DBAs 
(WD 0948+013, WD 2154-437) and two of our DBs 
(WD 0615-591, WD 0845-188) are in common proper 
motion or wide binar y systems (|Luvtenl 119491 : IWegner! 
[19731: lCaballero|[2009[ ). The lack of detectable Zeeman 
splitting implies that our targets also do not harbor sig- 
nifica nt (>100kG) magnetic fields (e.g., iKoester et al.1 
[20091) . 

The WDs in our sample have not been classified in the 
literature as potential members of the thick disk or halo 
stellar populations. We therefore make the assumption 
that all of these WDs belong to the thin disk - a neces- 
sary assumption for our analysis (see Section l3.3[) . Nearly 
all (>90%) of the SPY WDs that have been s tudied kine- 
matically are classified a s thin disk objects dPauli et al.l 
I21M iRlcliter et al.l[2007h . 



3. USING GRAVITATIONAL REDSHIFT TO DETERMINE A 
MEAN MASS 



We use the methods described in IFalcon et al.l (|2010[) . 
To summarize, we measure the apparent velocity of the 
Ha line for an ensemble of WDs. We correct these veloc- 
ities so that the WDs are in a comoving reference frame. 
The average apparent velocity then becomes the average 
gravitational redshift because random stellar radial ve- 
locities average out. By using the mass-radius relation 
from WD evolutionary models, we translate this average 
gravitational redshift to an average mass. 

In this work, we add two extensions to the method: 
(1) a different estimator of central location of the distri- 
bution better suited for small sample sizes (Section f3.5|) 
and (2) use of the He I 5876 A line for helium-atmosphere 
WDs (Section [3H). 



3.1. Gravitational Redshift 

In the weak-field limit, the general relativistic effect of 
gravitational redshift can be observed as a velocity shift 
in absorption lines and is expressed as 



cAA GM 



Be 
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where G is the gravitational constant, and c is the speed 
of light. In our case, M is the WD mass, and R is the 
WD radius. 

The apparent velocity of an absorption line is the sum 
of this gravitational redshift and the stellar radial veloc- 
ity: ^app = Vg + v T . These two components cannot be 
explicitly separated for individual WDs without an in- 
dependent v T measurement or mass determination. We 
break this degeneracy for a group of WDs by assuming 
that the sample is comoving and local. After we correct 
each i> app to the local standard of rest (LSR), only ran- 
dom stellar motions dominate the dynamics of our sam- 
ple. These average out, and the mean apparent velocity 
becomes the mean gravitational redshift: (v app ) = ( v g)- 

3.2. Velocity Measurements 

Collisional (Stark) broadening effects cause asymmetry 
in the wings of absorption lines for the hydrogen Balmcr 
series, making it difficult to measure a velocity c entroid 
(jShipman fc Mehanlll976t iGrabowski et al.lll987t) . How- 
ever, these effects are not significant in the sharp, non- 
LTE Ha and H0 Balmer li n e core s. We use both of these 
line cores in iFalcon et all ()2010[ ). In the SPY spectra 
of DBAs, though, we do not observe the non-LTE line 
cores, but the Ha line centers are still distinct. For Ha 
the pressure shift should remain very small (< 1 km s _1 ) 
w ithin a few A of the line center (jGrabowski et al1ll987f ). 

iReidl rfl996h measures the gravitational redshift of the 
DBA WD 0437+138 using the line center since he ob- 
serves no sharp line core. While cautioning that pressure 
effects may be signi ficant, he determines a high mass of 
0.748 ± 0.037 M Q . iBergeron et"aH (j2011[) determine a 
spectroscopic mass of 0.74 ± 0.06 Mq for this WD. For 
at least this one case, then, gravitational redshifts of the 
line center give a result consistent with the spectroscopic 
method. 

Van der Waals broadening due to neutral helium may 
also signficantly affec t line shapes belo w T e g =16,500 K 
(|Koester et alJl2005al iVoss et al.ll2007| ). We avoid this 
systematic by not including targets in that range of T a g . 

We measure i> app for each DBA in our sample by fitting 
a Gaussian profile to the Ha line center using GAUSS- 
FIT, a non-linear least-squares fitting routine in IDL. 
If multiple epochs of observation exist, we combine the 
measurements as a mean weighted according to the un- 
certainties returned by the fitting routine. Eight out of 
our 16 DBAs have two epochs of observation. Table [1] 
lists these measurements for the DBA sample including 
those for each epoch of observation and the final adopted 
values. 

By inspecting the v app measurements between epochs, 
we notice they differ by more than a typical measurement 
uncertainty; the mean difference between epochs of the 
same target is 9.20 km s^ 1 while the mean measurement 
uncertainty for all observations is 3.05 km s . Since we 
presume none of our targets to be in close binary sys- 
tems (Section 12. 1| . reflex orbital motion cannot be the 
culprit, and therefore it is evident that we are underes- 
timating our individual measurement uncertainties. To 
more accurately represent our adopted Svg. pp values for 
each target, we multiply all observation Sv app values by 
the factor we find above: the ratio of the mean difference 
in apparent velocity between epochs and the mean ap- 



parent velocity measurement uncertainty. For the DBA 
sample, this factor is 3.02. The final adopted values for 
<^ app listed in Column 3 in Table [T] include this adjust- 
ment, and the observation values in Column 8 are the 
values before the adjustment. 

For the DBAs and the DBs, we also measure i> app using 
the He I 5876 A line and list the results in Tables [2] and [3 
respectively. Here, too, we apply a multiplicative factor 
to all observation Sv app values to better represent our 
measurement uncertainties. These factors are 2.12 and 
2.49 for the two samples. The methodology for line center 
fitting is the same as used for the Ha Balmer line. 

3.3. Comoving Approximation 

We measure a mean gravitational redshift by assuming 
that our WDs are a comoving, local sample. With this 
assumption, only random stellar motions dominate the 
dynamics of our targets; this falls out when we average 
over the sample. 

For this assumption to be valid, at least as an approx- 
imation, our WDs must belong to the same kinematic 
population; in our work, this is the thin disk. We achieve 
a comoving group by correcting each measured w app to 
the k inematical LSR described by Standard Solar Mo- 
tion (jKerr fc Lvnden-Belll[l986l ). Column 4 in Tables Q] 
- [5] lists the LSR c orrections applied to each target. In 
IFalcon et al1(|2010t ). we show in detail that this is a suit- 
able choice of reference frame for a sample consisting of 
thin disk WDs. 

If the targets in our sample belong predominantly to 
the thick disk, for example, our chosen reference frame 
will fail to produce a group of objects that is at rest 
with respect to the LSR, introducing a directional bias 
into our measured (v app ). This is because the thick disk 
population lags behind the thin d isk as it rotates abo ut 
the Galactic center (~ 40 km s~ 1 : lGilmore et aLl [l989) 
In contrast to the work on DAs in lFalcon et al.l (|201C 
we lack the sample size to empirically demonstrate 
whether or not our WDs move with the LSR. Therefore, 
as we state in Section 12.11 we must assume our targets 
are thin disk objects. 

Our WDs must also reside at distances that are small 
compared to the size of the Galaxy so that system- 
atics introduced by the Galactic kinematic structure 
are not significant. Nine of our DBAs and twelve of 
our DBs have published dis tance determ i nations in the 
range of ~ 25 to ~ 260 pc (lRoutlvlll972t I Wegnen 119731: 
Koester et al.lll981t iFarihi et alJl2005b iCastanheira et al.l 



20061: iLimoges fc Bergeron! 120101 : IBergeron et all 120111 ) 



We assume the others are at di stances comparabl e to 
the rest of the SPY sample which iPauli et all (|2006ft de- 
termine spectroscopically to be mostly (>90 %) within 
200 pc with a mean of ~100pc. Over these distances, the 
velocity dispersion w ith varying scale height a bove the 
disk remains modest (jKuiiken fc Gilmorelll989|). and dif - 
ferential Galactic rotation is neglig ible dFich et alllT98l . 



3.4. Deriving a Mean Mass 

To translate the mean apparent velocity (u app ) to a 
mean mass, we invoke two constraints: (1) we need the 
mass-radius relation from an evolutionary model, and 
(2) since the WD radius does slightly contract during 
its cooling sequence, we need an estimate of the position 
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along this track for the average WD in our sample (i.e., 
a mean T c g). 

Our evolutionary models use Miie/M* = 10~ 2 for the 
helium surface-layer mass, and, since we are interested 
in WDs with helium-dominated atmospheres, we use 
Mh/Mj, = for t h e hyd rogen surface-layer mass. Sec 
I Montgomery et ahl (|1999() for a more complete descrip- 
tion of our models. Our dependency on evolutionary 
models is small. We are interested in the mass-radius re- 
lation from these models, and this is relatively straight- 
forward since WDs are mainly supported by electron de- 
generacy pressure, making the WD radius a weak func- 
tion of temperature. We estimate that varying the C/O 
ratio in the core affects the radius by less than 0.7%. 

3.5. Appyling the Method to Small Samples 

Besides atmospheric co mposition, th e diffe rence be- 
tween the WD sample in iFalcon et al.l (|2010f ) and the 
ones in this paper is sample size; the DA sample boasts 
449 targets while our DBA and DB samples contain 16 
and 20, respectively. For large numbers such as for the 
DAs, the (arithmetic) mean is the preferred estimator of 
central location of the distribution. For small numbers, 
however, this is not necessarily so. 

In lieu of the mean, we use the biweight central loca- 
tion estima tor and the biweigh t scale estimator as recom- 
mended by iBeers et al.l (|1990|) for a kinematic sample of 
our size. For a discussion on the biweight and its statis- 
tical properties of resistance, robustness, and efficiency 
as they rela t e to sa mple size and type of distribution, see 
IBeers et al.1 (fl990T ). 

The biweight central location estimator is defined as 



Cbi = M 



E w <i(i-«?) 



2\2 



(2) 



where Xi are the data, M is the sample median, and Ui 
are given by 

u.; = - . (6) 

cMAD v ' 

We set the "tuning constant" c = 7.0 (a sli ght increase 
from the recommended c = 6.0 listed in IBeers et al.l 
(|1990f l) so that no data rejection occurs. Notice that 
Cbi approaches the arithmetic mean in the limit where 
c — >• oo. MAD (median absolute deviation from the sam- 
ple median) is defined as 

MAD = median^ — M\). (4) 

We calculate Cbi iteratively, taking M as a first guess, 
substituting the calculated Cbi as the improved guess, 
and continuing to convergence. 
The biweight scale estimator is defined as 



£kui<i(s<-AQ a (i 



,2\4 



1/2 



£ w <i(i-^ 2 )(i-5^ 2 ) 



(5) 



Here we set c = 9.0 (within the e xpression for m) as 
recommended by IBeers et al.l (|1990f) . 

To determine the uncertainty of the biweight, we use 
a Monte Carlo approach. We adopt the standard devia- 
tion of biweight values from 30,000 simulated i; app distri- 
butions. Each simulated distribution randomly samples 
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Figure 1. Plot of apparent velocity u apP obtained from Hct versus 
v app obtained from the He I 5876 A line for the targets in our 
DBA sample. The intersection of (^app)Bl and (f a i p P )Bl (filled, 
red square) lies below the unity line; the two average values differ 
by 6.9 ±6.9 km s _1 . The dashed, bl ue line represents a theoretical 
bluesh ift of the He I line derived from Djmitrijcvic & Sahal- Brcc hot 
p990D for T = 16, 500 K and n e = 10 17 cm" 5 (see Section I3.6|l. 

from a convolution of two Gaussians. One represents the 
measured i; app distribution by using the biweight scale 
estimator Sbi (analogous to standard deviation) of the 
sample as the width of the Gaussian. The other accounts 
for the individual apparent velocity measurement uncer- 
tainties by drawing from a Gaussian with a width equal 
to each 5v ap p from the sample. 

From here on we use the subscript BI to denote the 
biweight central location estimator, Cbi, and biweight 
scale estimator, Sbi, as used in lieu of the arithmetic 
mean and standard deviation, respectively. 

3.6. Pressure Shifts of Helium Lines 

As described in Section [TJ helium line centroids are 
expected to shift due to pressure effects, and the mag- 
nitude and direction of the shift is thought to depend 
on the specific line. For that reason, we focus on mea- 
surements of a sing le helium line, He I 5876 A. Theo- 
retical calculations give a negative (blue) Stark shift of 
the He I 5876 A line for electron densities correspond- 
ing to WD photospheres and in the temperature range 
relevant to our investigation (e.g ., iGriem et al.l Il962t 
Dimitriievic fc Sahal-BrechotlrT990l ); the experiment by 
Heading et al.l (jl992P i confir ms the sign of the shift. 



We search for this expected blueshift using our sam- 
ple of DBAs, whose spectra show both the He I 5876 A 
line and the Ha Balmer line center that should not be 
affected by pressure shifts. Figure [T] shows a comparison 
of apparent velocities measured from these lines, and in- 
deed the intersection (filled, red square) of the average 
(biweight) values lies below the unity line. 

Using the calculated shifts from 

iDimitriievic fc Sahal-Brechotl (|1990f ). we plot a dashed 
line in Figure Q] indicating the deviation from unity 
for plasma conditions corresponding to temperature 



T = 16, 500 K and electron density 



10 17 cm" 3 . 



This should represent an overestimate of the expected 
shift based on theory. T e g = 16, 500 K is the lower 
limit for the WDs in our sample, and not only does the 
magnitude of the shift increase with decreasing temper- 
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ature, but it tracks exponentially so that the effect is 
more exaggerated at lower T. Furthermore, although 
n c = 10 17 cm -3 is typical for WD photospheres, we 
are concerned with absorption line centers. These are 
formed higher in the atmosphere where the density is 
lower and the shift should be less. As expected, the 
intersection of our measured biweight values lies between 
the unity line and the blueshift overestimate. 

The presence of hydrogen in a helium-dominated at- 
mosphere may also affect the atmospheric pressure and 
therefore electron density because of its relatively high 
opacity; in our temperature range, the opacity of hydro- 
gen can be an order of magnitu de greater than that of 
helium. From IVoss et al.l ([2007), the mean atmospheric 
hydrogen abundance for the WDs in our DBA sample is 
H/Hc= 10~ 329 with the highest abundance of any WD 
being H/Hc= 10~ 2 . At these abundances, this opacity 
effect should be negligible, and the dependency of the 
distribution of pressure shifts on electron density should 
only be due to the intrinsic mass distribution. 

The average apparent velocity of the DBAs in Figure [1] 
as measured from Ha is (v app )Bi = 40.8±4.7km s _1 ; the 
average velocity as measured from the He I 5876 A line is 
( u app)Bi = 34.0±5.0km s^ 1 . The measured difference of 
6.9±6.9km s _1 is blueshifted and exactly at the 1-sigma 
boundary of our measurement precision. 

With assumptions this measured blueshift can be used 
as a correction to apply to the average apparent ve- 
locity of DB WD samples in order to estimate an av- 
erage gravitational rcdshift. First, one must assume 
that DBAs arc not fundamentally different from DBs 
in a way that wo uld manifest as different mean masses. 
IVoss et all ()2007l ) detect various amounts of hydrogen 
in most (55%) of the helium-dominated WDs in their 
sample and find similar spectrosc opic mass distri b utions 
between the DBAs and DBs. iBergeron et al.l (|2011f ) 
also find no significant differences between the masses 
of the two groups. This provi des evidence to the idea of 
iWeidemann fe Koesterl (|1991[ ) that the DBA subclass is 
not distinct from its parent class but rather the obser- 
vationally detectable end of a continuous distribution of 
hydrogen abundances. 

Also, though our value is indeed a measure of 
the average blueshift due to pressure effects, this av- 
erage is not straight-forward. The magnitude and 
sign of the shift depend on t emperature and elec- 
tron density (or pressure) (e. g.. iKobilarov et al.l fl989t 
iDimitriievic fc Sahal-Brechotl 11990ft . and the atmo- 
spheres of our DBAs sample a distribution of these 
plasma conditions. For our measured value to be rep- 
resentative of the average blueshift for another sample 
of WDs, it is important that this sample have a similar 
T e g distribution. Assessing the similarity between sam- 
ples when applying this correction will deserve further 
scrutiny when a future analysis is performed which uses 
significantly larger sample sizes and hence higher preci- 
sion. 

We apply this correction so that we may, for the first 
time, estimate an average gravitational rcdshift for DB 
WDs in the field. The precision of this correction is given 
by the combined uncertainty of the method using the 
Ha Balmer line and of the one using the He line (i.e., 
^(^app)Bi and S(v^p P )bi added in quadrature; ~ 7km 
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Figure 2. Histograms of measured apparent velocities of Ha v app 
(dashed, filled) and of He I 5876 A ^ p e p (solid, unfilled) for 16 
DBAs. The bin size is 13 km s — . The dashed and solid curves are 
the Gaussian distributions used to determine Monte Carlo uncer- 
tainties for each sample. For H, (w a pp)Bl = 40.8 ± 4.7km s _1 and 
ctK pp )bi = !6-6km s" 1 . For He, (vg? p } B i = 34.0 ± 5.0km s^ 1 
and o"(^p C p)Bl = 18.6 km s _1 . 

s" 1 ) and can be improved by increasing the sample size 
of DBA WDs. 

4. RESULTS 
4.1. DBAs 
4.1.1. Analysis with the Ha Balmer Line 

We observe hydroge n absorption in 16 WDs classifed as 
helium-dominated bv IVoss et al.l (j2007l ). Figure [2] shows 
the distribution of u app measured from the Ha Balmer 
line. Table [1] lists these measured values, which are de- 
scribed in Section |3~21 

For this sample, (f app )Bi = 40.8 ±4.7 km s -1 . We find 
that this is different from (w app ) of DAs (32.57km s _1 ) 
at the 91% confidence level. 

Given the small sample size, it is possible that a sys- 
tematic uncertainty is skewing our result. For example, 
by chance kinematic substructures may not average out 
with only 16 targets as they would for a larger sample. 
As a check, we perform the following test: let us assume 
that the SPY DBAs have the same intrinsic mean mass as 
the SPY DAs. If wc randomly pick 16 appar ent velocities 
from the measured v app distribution from iFalcon et al.l 
(|2010f ). which contains 449 objects, how often is the 
(^app)Bi of these 16 greater than (i> app )Bi = 40.8 km s _1 , 
our result for DBAs? Performing the random selection 
30,000 times, we find that this false detection occurs 9% 
of the time. The synthetic biweights reproduce the (^ app ) 
of DAs with a standard deviation of 4.17km s _1 . There 
is a rare but not insignificant possibility that our small 
sample size is fooling us. 

Using s pectro scopically determined T e g from 
IVoss et all (l200l . (T eff ) B i = 17,890 ± 250 K. We 
use the mass-radius relation from our evolutionary 
models and interpolate to find (M)bi = 0.711q'q5 M©. 
This is larger than the mean mass for normal DAs de- 
termined using the same gravitational redshift method 
(0.647 j^g" M Q ; IFalcon et al.llMol) . 

Though suggestive of a real difference in mean mass 
between the samples, this is not a stringent result. It 
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remains plausible that no intrinsic difference exists. We 
can say, however, that the difference between (u a p P ) B i of 
the DBAs and (w a pp) of normal DAs is no larger than 
9.2 km s _1 , at the 95% confidence level. With a note of 
caution, because of the sublctics in translating a mean 
apparent velocity to a mean mass, this corresponds to 
roughly 0.10 M Q . 

Our (M)bi for DBAs is also larger than 0.62±0.02 M , 
the biweight of the spectroscopic mas s determinations for 
these targets from IVoss et aLl (|2007t ) . The value for the 
mean is the same as the biweight. We estimate the un- 
certainty u sing the same met hod as for our work, except 
that, since IVoss et ail (120071 ) do not list individual un- 
certainties, we assign each mass an uncertainty equal to 
the standard deviation of the mass distribution. The bi- 
weight (and mean) spectr oscopic mass of the 22 DBAs 
with > 16, 500 K from iBergeron et alJ (|2011[ ) agrees 
well at 0.67 ± 0.02 M Pl . For this value w e use the mass 
uncertainties from IBergeron et aLl (|201lD . 

4.1.2. Analysis with the He I 5876 A Line 

By measuring v^ pp - instead of that from the Ha 
line - for the targets in our DBA sample, we find 
( u a i p°p)Bi = 34.0 ± 5.0km s" 1 . Table M lists individ- 
ual measurements. As discussed in Section 13.61 this 
differs from the Ha result by 6.9 ± 6.9 km s — 1 . It 
lies between t he unity line and the theoretic a l blue shift 
derived from iDimitriievic fe Sahal-Brechotl (|1990f ) for 
plasma conditions corresponding to T = 16, 500 K and 
n c = 10 17 cm -3 . This temperature is the lower limit of 
our sample, and this electron density is typical of WD 
photospheres. 

4.2. DBs 

We now use the He I 5876 A line to perform the grav- 
itational rcdshift method on DBs. We measure uH?L for 

app 

20 such WDs and find (v^ p ) B i = 36.0 ± 5.6 km s _1 . 
Table [3] lists individual measurements. After applying 
the correction for the blueshift due to pressure effects 
measured in Section E21 (Vapp)Bi = 42.9 ± 8.9km s -1 . 
This uncertainty comes from adding the measured uncer- 
tainty and that of the co rrection in quad rature. Using 
(T eff ) B i = 20, 570 ± 660 K (jVoss et al.ll2007h with the cor- 
rected (Wapp)Bi, we determine (M) B i = O.74t°;ogM . 

This is different from (v a p P ) of DAs at the 75% con- 
fidence level. The difference between {v^ )bi (with the 
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blueshift correction) of the DBs and 



,) of DAs is 



< 11.5 km s 1 , at the 95% confidence level. This corre- 
sponds to roughly . 12 Mp . 

From IVoss et alJ EMM) , we find (M*) B i = 0.58 ± 
0.02 Mq for these targets (M* denotes spectroscopic 
mass); the mean is 0.59 ± 0.02 Mq. The biweight spec- 
troscopic mass o f the 35 DBs with T e g > 16, 500 K from 
IBergeron efaTl (|2011l) is 0.63 ± 0.01 Mq] the mean is 
0.64 ± 0.01 M Q . The mean of the 82 DBs with T eS > 
16, 000K from SP SS and whose spect ra have S/N> 20 is 
0.646 ± 0.006 Mq (jKepler et al.ll20lTit) . 

4.3. DBAs and DBs 

As mentioned in Section 13.61 let us assume there are 
no fundamental differences between DBAs and DBs that 



CO 
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-20 20 40 60 BO 100 
H Apparent Velocity (km s~') 

Figure 3. Histogram of measured apparent velocities of the He I 
5876 A line v^ p for 36 DBAs and DBs (16 and 20, respectively) 
with a bin size is 15 km s — 1 . The dashed curve is the Gaus- 
sian distribution used to determine Monte Carlo uncertainties. 
( 1 'app)Bl = 36.5 ± 3.9km s~ 1 (before applying the correction dis- 
cussed in Section 13.611 and ct(ii^ p c p )bi = 22.4km s — 1 . 

would manifest as different mean masses. We can then 
apply the correction for the blueshift due to pressure ef- 
fects, and we can improve both the precision and accu- 
racy of our mean value determinations by combining the 
two samples. We do not discard the possibility, however, 
that the two groups are indeed fundamentally different. 

Using the apparent velocity of the He I 5876 A line 
for the combined sample of DBAs and DBs, we find 
( w app)Bi = 36.5 ± 3.9 km s -1 . Figure [3] shows the dis- 
tribution of v^ pp . Applying the blueshift correction, 
^a P p>Bi = 43.4 ± 7.9km s" 1 . Using (T ofr ) BI = 19, 260 ± 



390 K, we determine (M) B i = 0.74+°;°^ M Q . 

This is different from the mean w app of DAs at the 
82% confidence level. The difference between {w a i p C p ) B i 
of this combined sample of DBAs + DBs and (i) ap p) of 
DAs is < 7.7 km s _1 , at the 95% confidence level. This 
corresponds to roughly 0.08 Mq. 

(M*) B i = 0.60 ± 0.02 Mm for thes e targ ets using the 
spectroscopic masses from IVoss et alJ (|2007f ). This is sig- 
nificantly lower than our determination. The biweight 
(and mean) spectroscopic mass of the 57 DBAs and 
DBs with T eS > 16, 500 K from IBergeron et all (p?5Tlh 
is 0.65 ± 0.01 Mq. 

We list our derived (u app ) B i in Tabled denoting which 
samples use the He I 5876 A line in lieu of the Ha Balmcr 
line. Since the (^app^i values are entirely model inde- 
pendent, we list them apart from the (Af) B i we deter- 
mine. These are in Table [5] For co mparison, we list th e 
corresponding results for DAs from I Falcon et alJ (|2010f ) . 

5. CONCLUSIONS 
We measure the apparent velocity (u app ) of the He I 
5876 A line for a sample of DBAs and compare it to 
that of the Ha Balmer line. We find a difference 
of 6.9 ± 6.9 km s _1 in the average apparent velocities 
from the two line species, which we attribute to the 
blueshift of this He line due to p ressure effects (e.g., 
IDimitriievic fe Sahal-Brechotl I1990D averaged over the 
sample. With assumptions one can apply this measured 
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blueshift as a correction to other samples of helium- 
atmosphere WDs. We do so in order to investigate the 
average gravitational redshift of a sample of DBs for the 
first time. 

Following the g ravitational redshift method from 
iFalcon et all (|2010l ) . but using biweight estimators which 
are better suited for small sample sizes, we find (v s )bi = 
(u app ) B i = 40.8 ± 4.7 km s" 1 for 16 DBAs with T cff > 
16, 500 K from SPY. We translate this (w ap p}Bi to a 
mass: (M) B i = 0.71^;^ M . Though different from 
the (w a pp) of DAs, 32.57 km s -1 , at the 91% confidence 
level and suggestive of a larger DBA mean mass than DA, 
OM7±o° {lM Q (jFalcon et al.ll2010h . this is not a strin- 
gent result. We emphasize that the difference between 
(fapp)Bi of the DBAs and (w a pp) of normal DAs is no 
larger than 9.2 km s _1 , at the 95% confidence level; this 
corresponds to roughly 0.10 Mq. Our (M)bi for DBAs 
is also larger than the average of the s pectroscopic mass 
determinations for these targets from IVoss et al.1 ()2007l ) 
at 0.62 ± 0.02 Mq. It agrees with the average mass of 
the 22 DBAs with T off > 16, 500 K from iBergeron et al.l 
d20ll at 0.67 ± 0.02 M . 

We use the He I 5876 A line to conduct the first gravita- 
tional redshift investigation of a group of WDs without 
visible hydrogen lines. For 20 DBs from SPY, we find 
( u app)Bi = 42.9 ± 8.9 km s _1 after applying the correc- 
tion for our measured blueshift due to pressure effects. 
We determine (M) B i = 0.74+^f M . The difference be- 
tween (Wapp)Bi of the DBs and (w ap p} of DAs is < 11.5 km 
s _1 , at the 95% confidence level; this corresponds to 
roughly 0.12 Mq. The (M)bi is much larger than the 
average of the spectroscopic mass determinations from 
SPY, 0.58 ± 0.02 M Q . It is slightly larger than the aver- 
age spectroscop ic mass of DBs with T e g > 16, 500 K from 
IBergeron etafl (|2011l ). 0.63 ± 0.01 Mq, and the mean of 
DBs with T fi ff > 16 000 K from SDSS, 0.646 ± 0.006 M 
(Kepler et al.ll20Toh . 

Combining our DBA and DB samples to group all WDs 
with helium-dominated atmospheres, we find (v^ p c p )bi = 
43.4 ± 7.9 km s _1 (after the correction) and determine 
(M)bi = O.74t°;^M . This differs from the (v app ) 
of DAs at the 82% confidence level. The difference 
between (Vapp) BI °^ * ne DBAs + DBs and (i> app ) of 
DAs is < 7.7km s" 1 (roughly 0.08M o ), at the 95% 
confidence level. Our (M)-qi is much larger than the 
average spectroscopic mass of these targets from SPY 
at 0.60 ± 0.02 Mq and slightly larger than the average 
spectroscopic mass of the DBAs an d DBs with T e g > 
16,500K from IBergeron et al.l poll . 0.65 ± 0.01 Mq. 

Given the small sample sizes, it is possible that system 
uncertainties are skewing our results due to the potential 
of kinematic substructures that may not average out. We 
estimate this to be unlikely, but a larger sample size is 
necessary to rule out these systematics. 
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Table 1 

H Apparent Velocity Measurements for DBAs 





Adopted 


Date 


Time 


LSR 


Observation 




^app 


dfapp 








Correction 


^app 


dfapp 




(km s J 


(km s J 


CUT) 


l u 1 ) 


(km s ) 


(km s 1 ) 


(km s — 


HE 0025-0317 


36.76 


5.45 


2000 


.07.17 


07:52:34 


27.319 


36.76 


1.80 


HE 0110-5630 


42.04 


8.01 


2002 


,09.25 


07:28:16 


-16.836 


42.04 


2.65 


WD 0125-236 


55.57 


13.45 


2002 


,09.30 


04:23:09 


-6.499 


55.57 


4.45 


WD 0921+091 


46.71 


7.90 


2001 


,04.07 


00:55:54 


-35.250 


53.62 


2.49 








2002 


12.28 


07:48:25 


10.131 


42.19 


2.01 


WD 0948+013 


30.05 


6.12 


2001 


,01.10 


01:31:21 


-34.417 


38.21 


4.20 








2003 


,01.17 


04:03:53 


4.365 


27.75 


2.23 


WD 1149-133 


48.92 


3.98 


2000 


,07.13 


23:40:32 


-34.214 


53.48 


3.71 








2000 


,07.16 


23:55:06 


-33.761 


47.18 


2.29 


HE 1207-2349 


31.14 


4.02 


2002 


,02.23 


07:41:20 


12.910 


43.82 


6.79 








2002 


,02.24 


07:11:53 


12.571 


30.51 


1.52 


EC 12438-1346 


13.50 


16.32 


2000 


,07.15 


00:27:35 


-31.791 


13.50 


5.80 


WD 1311+129 
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5.28 


2001 


,06.20 


01:47:22 


-20.994 
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2.87 
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,01.18 


08:36:37 


33.896 


6.55 


4.60 


HE 1349-2305 
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5.91 


2000 


,07.15 


02:59:13 


-28.698 


15.58 


1.95 


WD 1421-011 


56.32 


11.11 


2001 


,08.16 


00:32:05 


-19.271 


56.32 


3.67 


WD 1557+192 


52.12 


4.68 


2002 


,04.23 


06:37:23 


22.663 
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1.55 


WD 1709+230 


72.70 


1.68 


2002 


,09.04 


00:10:33 


-1.552 


71.41 


3.30 








2002 


,09.21 


00:32:10 


-1.091 


73.80 


3.05 


WD 2154-437 


41.92 


4.20 


2000 


,06.04 


06:01:21 


25.025 


41.92 


1.39 


WD 2253-062 


41.47 


4.32 


2000 


,06.01 


08:50:29 


37.290 


37.64 


3.38 
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,06.06 


08:07:23 
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44.02 
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HE 2334-4127 
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39.11 


2.67 
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Table 2 

He Apparent Velocity Measurements for DBAs 



Target 


Adopted 


Date 


Time 


LSR 


Observation 




? .-tie 

app 

(km s 


™app 

) (km s _1 ) 






Correction 


„Hc 
app 

(km s _1 ) 


du app 

(km s _1 ) 




(UT) 


(UT) 


(km s" 1 ) 


HE 0025-0317 


32.90 


9.63 


2000.07.17 


07:52:34 


27.319 


32.90 


4.54 


HE 0110-5630 


39.01 


7.55 


2002.09.25 


07:28:16 


-16.836 


39.01 


3.56 


WD 0125-236 


42.83 


4.43 


2002.09.30 


04:23:09 


-6.499 


42.83 


2.09 


WD 0921+091 


37.67 


4.63 


2001.04.07 


00:55:54 


-35.250 


40.38 


2.55 








2002.12.28 


07:48:25 


10.131 


33.71 


3.08 


WD 0948+013 


10.98 


1.33 


2001.01.10 


01:31:21 


-34.417 


12.76 


2.76 








2003.01.17 


04:03:53 


4.365 


10.48 


1.46 


WD 1149-133 


42.45 


0.18 


2000.07.13 


23:40:32 


-34.214 


42.13 


6.68 








2000.07.16 


23:55:06 


-33.761 


42.51 


2.68 


HE 1207-2349 


26.12 


8.75 


2002.02.23 


07:41:20 


12.910 


20.72 


3.80 








2002.02.24 


07:11:53 


12.571 


33.22 


4.36 


EC 12438-1346 


15.56 


3.04 


2000.07.15 


00:27:35 


-31.791 


15.56 


1.43 


WD 1311-1-129 

v v \—j tun i J 


22.72 


2.69 


2001.06.20 


01:47:22 


-20.994 


21.59 


2.15 








2003.01.18 


08:36:37 


33.896 


25.91 


3.61 


HE 1349-2305 


10.49 


3.72 


2000.07.15 


02:59:13 


-28.698 


10.49 


1.75 


WD 1421-011 


66.97 


8.23 


2001.08.16 


00:32:05 


-19.271 


66.97 


3.88 


WD 1557+192 


61.57 


3.28 


2002.04.23 


06:37:23 


22.663 


61.57 


1.55 


WD 1709+230 


ni a a 
yl.4D 


3.78 


2002.09.04 


00:10:33 


-1.552 


on ni 
oy.y 1 


l.oo 








2002.09.21 


00:32:10 


-1.091 


96.09 


2.73 


WD 2154-437 


30.39 


4.40 


2000.06.04 


06:01:21 


25.025 


30.39 


2.07 


WFl 99^ 0fi9 


39.21 


3.64 


2000.06.01 


08:50:29 


37.290 


37.50 


1.36 








2000.06.06 


08:07:23 


37.293 


43.11 


2.06 


HE 2334-4127 


40.83 


6.33 


2002.09.14 


02:03:49 


-9.502 


37.32 


2.17 








2002.09.15 


01:27:02 


-9.859 


46.56 


2.77 








Table 3 














He Apparent Velocity Measurements for DBs 






Target 


Adopted 


Date 


Time 


LSR 


Observation 




v tlc fotlc 
app app 

(km s _1 ) (km s _1 ) 






Correction 


app 

(km s~ 1 ) 


<5u Hc 

app 

(km s- 1 ) 




(UT) 


(UT) 


(km s" 1 ) 


MCT 0149-2518 


62.43 


7.83 


2002.09.27 


08:32:10 


-3.801 


55.45 


2.56 








2002.09.30 


04:09:09 


-4.575 


66.83 


2.03 


WD 0158-160 


47.19 


3.20 


2002.09.19 


03:40:58 


0.897 


48.67 


1.61 








2002.09.20 


03:55:16 


0.462 


43.71 


2.46 


WD 0249-052 


46.41 


2.16 


2002.09.20 


07:09:51 


9.064 


50.32 


3.97 








2002.12.31 


02:03:00 


-33.912 


45.81 


1.55 


HE 0308-5635 


63.79 


6.87 


2002.09.14 


05:54:31 


-12.226 


63.79 


2.75 


WD 0349+015 


8.96 


2.21 


2003.01.24 


03:40:30 


-38.561 


8.96 


0.88 


HE 0417-5357 


2.19 


8.48 


2000.07.15 


10:16:58 


-8.104 


-6.21 


3.25 








2000.07.18 


09:37:59 


-8.014 


6.47 


2.31 


HE 0420-4748 


57.29 


0.06 


2000.07.15 


10:30:58 


-6.686 


57.37 


2.83 








2000.07.18 


10:19:41 


-6.498 


57.26 


1.52 


HE 0423-1434 


55.50 


5.98 


2002.09.23 


06:50:16 


4.515 


45.45 


4.07 








2002.12.07 


08:17:12 


-22.990 


57.28 


1.71 


WD 0615-591 


1.94 


0.61 


2000.09.14 


09:26:25 


-13.891 


1.322 


2.83 








2001.04.08 


01:04:56 


-21.666 


2.251 


2.00 


WD 0845-188 


72.10 


5.00 


2001.12.18 


07:32:37 


2.566 


70.29 


1.79 








2001.12.29 


07:17:47 


-0.614 


79.02 


3.50 


WD 1252-289 


16.99 


0.59 


2000.04.21 


04:49:17 


-8.262 


16.83 


1.46 








2000.05.19 


01:56:15 


-20.292 


18.10 


3.85 


WD 1326-037 


51.82 


6.61 


2001.01.14 


09:11:41 


33.888 


59.85 


1.91 








2001.02.02 


09:27:35 


31.681 


49.09 


1.11 


WD 1428-125 


48.55 


3.51 


2000.07.06 


02:51:29 


-20.736 


53.20 


5.01 








2000.07.13 


02:26:45 


-22.017 


47.22 


2.68 


WD 1445+152 


-4.35 


3.48 


2002.04.23 


05:32:37 


11.695 


-4.35 


1.39 


WD 1542+182 


27.89 


4.49 


2002.04.23 


08:22:15 


20.960 


27.89 


1.80 


WD 1612-111 


45.50 


2.48 


2000.06.06 


06:35:11 


6.461 


47.37 


2.68 








2000.06.08 


02:30:50 


5.989 


43.85 


2.52 


WD 2144-079 


9.97 


3.90 


2000.05.17 


09:50:49 


40.410 


11.93 


1.75 








2000.05.19 


09:58:26 


40.360 


6.080 


2.47 


WD 2234+064 


19.21 


4.19 


2001.06.18 


09:59:17 


37.415 


14.75 


4.18 








2001.09.01 


01:22:15 


12.083 


17.27 


2.54 








2002.09.24 


04:30:27 


0.547 


23.40 


2.75 


WD 2354+159 


44.02 


1.21 


2002.09.04 


06:44:23 


18.606 


43.46 


1.43 








2002.09.13 


04:25:56 


14.692 


45.33 


2.19 


WD 2354-305 


34.30 


10.00 


2000.07.15 


05:28:58 


20.072 


34.30 


4.01 



Falcon, Wingct, Montgomery, & Williams 



Table 4 

Biwcight Apparent Velocities 



Sample 


Number of WDs 


(^app)BI 

(km s- 1 ) 


<5(fapp)BI 

(km s _1 ) 


0-(Uapp)BI 

(km s- 1 ) 


(Sv ap p) 
(km s- 1 ) 


(M/R) m 

(Mq/Rq) 


S{M/R) BI 
(Mq/Rq) 


DBA 


16 


40.8 


4.7 


16.6 


7.1 


64.2 


7.4 


DBA (He) a 


16 


34.0 


5.0 


18.6 


4.7 


53.4 


7.8 


DB (He) b 


20 


42.9 


8.9 


23.9 


1.1 


67.4 


13.9 


DBA+DB (He) b 


36 


43.4 


7.9 


22.4 


1.5 


68.9 


12.5 


DA C 


449 


32.57 


1.17 


24.84 


1.51 


51.19 


1.84 



a Sample before pressure shift correction described in Section 13.61 
b Sample with applied pressure shi ft correction. 

c Main sample of normal DAs from(Jalcon ct al. (201(J). The values listed in this row contain arithmetic means, 
standard deviation, and uncertainties and no biweights. 



Table 5 

Biwcight Masses 



Sample 


Number of WDs 


(^app)BI 

(km s- 1 ) 


S (^app)BI 

(km s — 1 ) 


(T c ff )bi 
(K) 


o"(T cf r) B i 
(K) 


(AObi 
(Mq) 


6{M) m 
(Mq) 


DBA 
DB (He) 
DBA+DB (He) 
DA a 


16 
20 
36 
449 


40.8 
42.9 
43.4 
32.57 


4.7 
8.9 
7.9 
1.17 


17890 
20570 
19260 
19400 


1040 
3280 
2770 
9950 


0.71 
0.74 
0.74 
0.647 


+0.04 
-0.05 
+0.08 
-0.09 
+0.07 
-0.08 
+ 0.013 
-0.014 



a Main sample of normal DAs from Falcon ct al. ( 2010). The values listed in this row contain arithmetic 



means, standard deviation, and uncertainties and no biweights. 



